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Abstract 

In this paper we construct a sequence of eigenfunctions of the "quantum Arnold's 
cat map" that, in the semiclassical hmit, show a strong scarring phenomenon on the 
periodic orbits of the dynamics. More precisely, those states have a semiclassical 
limit measure that is the sum of 1/2 the normalized Lebesgue measure on the torus 
plus 1/2 the normalized Dirac measure concentrated on any a priori given periodic 
orbit of the dynamics. It is known (the Schnirelman theorem) that "most" sequences 
of eigenfunctions equidistribute on the torus. The sequences we construct therefore 
provide an example of an exception to this general rule. Our method of construction 
and proof exploits the existence of special values of h for which the quantum period 
of the map is relatively "short", and a sharp control on the evolution of coherent 
states up to this time scale. We also provide a pointwise description of these states 
in phase space, which uncovers their "hyperbolic" structure in the vicinity of the 
fixed points and yields more precise localization estimates. 
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1 Introduction 



One of the main problems in quantum chaos is the understanding of the semiclassical 
behaviour of the eigenfunctions of quantum dynamical systems having a chaotic classical 
limit. The main theorem in this context is the Schnirelman theorem |5cl ICdVl \Z1[ iHMRt 
IBouDBj . It roughly states that "most" eigenfunctions equidistribute on the available phase 
space in the classical limit. This leaves open the question of the existence of exceptional 
sequences of eigenfunctions with a different limit. In the case of "hard chaos" (uniformly 
hyperbolic systems), numerical computations have shown the presence of "scars" on certain 
eigenfunctions jHej . i.e. a visual enhancement of the wavefunction on an unstable periodic 
orbit. Up to now all theories of this phenomenon have required some kind of averaging 
over a (semiclassically large) set of eigenfunctions |Bog IBer| IHel IKH] . In addition, scarring 



is often described in the physics literature as a weak type of localization, compatible with 
Schnirelman's (measure-theoretic) equidistribution, as opposed to "strong scarring" |R Sj . 
which implies that the limiting measure has a component supported on a periodic orbit 
and therefore does not equidistribute. We show in this paper that, for the quantized 
"Arnold's cat map" , strongly scarred sequences do indeed exist for any periodic orbit (more 
generally, for any finite union of periodic orbits). This is, to the best of our knowledge, the 
first example of this kind in hyperbolic systems. A construction of exceptional sequences of 
eigenfunctions not equidistributing in the semiclassical limit was recently announced |CKSj 
for the quantization of certain ergodic piecewise afiine transformations on the torus, but 
these do not correspond to "scars" since the systems in question have no periodic orbits. 

Our construction is based on intuitively clear ideas that we now briefly sketch. For 
unfamiliar notation, we refer to Sections |2JH Precise statements of our results will be 
given below. 

Let M G SL(2,Z) be a hyperbolic automorphism of the 2-dimensional torus T and 
M its quantization on the A^- dimensional quantum Hilbert space Ti-N^g, where 211 hN = 1. 
We will construct strongly scarred quasimodes of M that, for certain values of A^, will be 
shown to be eigenfunctions. For that purpose we will use three ingredients. First, the 
time-energy uncertainty relation in the following simple form (T G N, G M): 

T-l 

II (M - e'^I) e"''^*M* || = || e-^'^^^M^^ - I ll< 2. (1) 

t=-T 

Second, precise estimates on intuitively clear phase space localization properties of coherent 
states. Third, a remark on the quantum period of M |BonDBl| (Section |HI). 

Let xo,xi = Mxo, . . . ,Xr = M'^xq = xq be a periodic orbit of period r of M. Let 
|xo,co,6') be a "squeezed" coherent state in TiNfi centered on the point xq and consider 
M*|a;o, Co, 6) for t E Z. Note first that this state is still a squeezed coherent state and that, 
for small enough t, it is localized around Xt- In fact, the support of the Husimi function of 
this state is an ellipse stretched along the unstable direction of the dynamics through the 
point Xt, with its major axis roughly of size \/^e^*, where A is the (positive) Lyapounov 
exponent of the dynamics (Section 0]). Introducing the Ehrenfest time T = the 
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support is therefore microscopic as long as t < (1 — e)T/2. For longer times, between T/2 
and T, the support of the Husimi function of M*|xo, Cq, 9) starts to wrap around the torus 
and it was shown in |BonDBl| that it equidistributes on that time scale. 
We shall consider the "discrete time quasimode" 

T-l 4 

= E e-^^*M*|xo, 60, 0) = Y^ (2) 

t=-T j=l 

and its "components" 

-T+if-l 

lO= E e-^*M*|xo,co,^). (3) 



t=-T+(j-l)f 



We note that similar states were considered before in the study of scars, see for instance 
|dPBB[ iKHj and references therein. We shall introduce a "continuous time" version 1$^°"^*) 

-T -T/2 T/2 T , 



l^l> l^2> l'^3> l'^4> 



' erg ' loc ' erg 



!<!>> 

Figure 1: Partition of the time interval [— T, T] into four equal parts, and of the quasimode 
into corresponding components. 



of those quasimodes later. We will write | $0) in statements true for both the discrete and 
continuous time quasimodes. 

Let us for simplicity concentrate on the case where = 0, r = 1. Our crucial technical 
estimate (Section EJ-Proposition H)) says that there exists C > so that 

(co,^ I M* I 5o,^) = -^=L== + m, with \I{t)\ < Ce-"(T-^). (4) 

'cosh(At) 



This implies rather easily (Proposition |21) the existence of a smooth, strictly positive func- 
tion 5*1 (0, A) so that 

(<l>^ I <l>^) ~25i(0,A)T. 
Using (dJ one concludes readily that 



(M-e'^i)|<|.^)J|< J— -^(1 + -^^) , (5) 



81(0, A)T V Si(0,A)T 
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justifying the name "quasimode" . Here we used the notation \ip)n = \/ {'>p\'4') for any 
non-zero G 7iM,e- 

To analyze the phase space properties of the above quasimodes, we first show as a 
further consequence of (jH) that the four states |$j,<^) have the same norm, asymptotically 
proportional to \/T as h goes to and that they are asymptotically orthogonal in the 
semi classical limit. In fact, this is easily understood intuitively by noting for example that 
the Husimi function of is supported along the stable manifold of the periodic orbit, 

and that of |$4,</)) along the unstable one, so that they have essentially disjoint supports, 
which is at the origin of their orthogonality. To put it differently, since the unstable and 
stable manifolds intersect at homoclinic points, our results show that the contribution of 
these intersections in the phase space integral expressing the overlap ($i,<^|$4,0) is small for 
small h. Note that although the homoclinic interferences do not contribute significantly to 
the above integral, they are nevertheless clearly visible on the pointwise behaviour of the 
Husimi distribution of which is represented in Figure|21and that will be further studied 
in Section El (for "continuous time" quasimodes). The pointwise estimates obtained there 
will show that the Husimi density concentrates along "classical hyperbolas" asymptotic to 
the stable and unstable manifolds; they will at the same time provide estimates on the 
rate of convergence to the limit measure, as well as other localization indicators (namely, 

norms of the Husimi density). 




Figure 2: Husimi distribution of the state |$0)„, constructed for the cat map (PT|) on 
the orbit of period 3 starting from xq = (0,0.5). The quantum parameters read = 
l/(27r^) = 500, = 0. (a): 3D plot on a linear scale, (b): 2D plot in logarithmic scale 
(darker = higher values). 

It is furthermore clear from the previous discussion on the phase space localization 
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properties of the evolved coherent states that |$i,<^) and |$4,</,) are sums of states that 
equidistribute on the torus, whereas |$2.(/.) and {^3^ are sums of states that locahze on 
the periodic orbit. One therefore expects (and we shall prove in Sections EHZj) that 



and that 



lim n{^j,^\f\^j,<f>)r 



lim „($j,0|/|$,,<^)„ 



f{x)dx ifj = l,4. 



1 ^"^ 

i=0 



if J = 2, 3. 



Here / is either the Weyl or anti-Wick quantization of / G C°°(T). In other words, the 
Wigner and hence also the Husimi function of \^2,3,<t>) converge (weakly) to the Dirac 
measure on the periodic orbit, whereas the ones of |$i,4,0) equidistribute, i.e. converge to 
the Lebesgue measure. This suggests grouping these states two by two, defining: 



6) = + |$4,<^) and |$ioc,</,) = 1^2,0) + |$3,0)- 



(6) 



Using the above information we shall finally prove (Propositions [7| and IT^ that, for 
any G [-7r,7r], 



lim „($^|/|<l><^). 



f{x)dx + 



J2 



j=0 



(7) 



In other words, the semiclassical limit measure of the sequence of quasimodes |$(/,)r 
measure 



is the 



1^ 1 
-ax H — 
2 2 



1 

r 



This shows that the quasimodes are strongly scarred. 

We then conclude using a particular property of the quantum period of M. We recall 
that the quantum cat map M has an h dependent "quantum period" P, i. e. = e~"^I 
for some (f G [0, 27t[. The eigenvalues of M on Ti.N,9 are therefore all of the form e~"'^j, with 
(pj = ip/P + 27rj/P, j = 1, . . . , P. Note that P plays the role here of the Heisenberg time 
of the system, since Acpj ~ 1/P. Since, for general h, the quantum period P is of order 
|Kej . it is considerably longer than the Ehrenfest time T, which grows only logarithmically 
in h~^. Nevertheless, developing an argument in |BonDBl] . we will show that, for any 
hyperbolic matrix in SL(2, Z) there exists a subsequence {hk)ken of values of h tending 
to zero for which P = 2T + 0{1) (see also |KR2j ). For those values the Heisenberg and 
Ehrenfest times of the system coincide and the |$</,)n therefore constitute a sequence of 
eigenfunctions of M that strongly scar, provided = (pj for some j G {1 . . . P}. It should 
be noted that, for the values of h considered, the number of distinct eigenvalues (pj is of 
order | lnh\, so that the eigenvalue degeneracy is very large, namely of order {h\ ln^|)~^. 

Our main result can finally be summarized as follows: 
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Theorem 1. Let M and {hk)keN be as above. Let < /? < 1/2 and let V = {xq, . . . ,Xr-i} 
be a periodic orbit of M. Then there exists a sequence {'4'jk)ke'M of eigenfunctions of M on 
'H.Nkfi with the property that, for all f G C°°(T^), 

lim ni^klfmn = /3-y\ fi^j) + (1 - /9) /" fix) dx. (8) 

Our result helps to complete the picture of the semiclassical eigenfunction behaviour of 
quantized toral automorphisms known to date. Indeed, beyond the general Schnirelman 
theorem for these models [BouDBj the following results are known. First, suppose M 
is of "checkerboard form", meaning AB = = CD mod 2. Then all eigenfunctions of M 
semiclassically equidistribute, provided one takes the limit along a density one subsequence 
of values of N |KR2j . for which the quantum period is larger than Note that this 

sequence excludes the values Nk for which the period is very short. Second, it is shown 
in |KRH IMej that for such M there exists a basis of eigenfunctions that equidistribute 
as tends to infinity, without restrictions on A^. This basis is constructed as a common 
eigenbasis for M and its "quantum symmetries" , which are shown in |KRlj to be sufficiently 
numerous to drastically reduce (if not to lift) the degeneracies of the eigenvalues. Finally, 
one may wonder if it would be possible to construct a sequence of eigenfunctions of M that 
has as a limit measure 

1 ""^ 

with /3 > 1/2. It is proven in |FNlj that this is impossible, so that the above quasimodes 
are in a sense maximally localized (the bound /? > {Vb — l)/2 = 0.62 had been previously 
obtained by jBonDB2j). 



2 Linear dynamics on the plane 

In this section we recall some known results we will need in the sequel. For details not 
given here we refer to [Fj. 



2.1 Classical linear flow 

The most general quadratic Hamiltonian on is (a,/5,7 G R): 

H{q, p) = ^aq^ + + iqp- (9) 

Assuming 7^ > a/?, H generates a hyperbolic flow x{t) = {q(t),p(t)) on R^, given by 
x{t) = M{t)x{0) (t G R), where for each t ^ 0, M(t) is a hyperbolic matrix in SL(2,R). 
Explicitly, for t = 1 

M^^M{1)=( ^ ^)gSL(2,R), (10) 
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i.e. AD-BC = 1, and 

A = cosh A + ^ sinh A _B = ysinhA 
C = — J sinh A D = cosh A — ^ sinh A 

where A = — ajS > is the Lyapounov exponent. Note that M has two real eigenval- 
ues e^'^ and hence two real eigenvectors corresponding to an unstable and a stable direction 
for the dynamics. They have respective slopes s+ = tan S-_ = tan Clearly, any hy- 
perbolic matrix M G SL(2, M) with TrM > 2 is of the above form for a unique a, (3, 7 (the 
case TrM < —2 is treated by using the map — M). The expressions in fllOj) - (fTT|) still make 
sense in the elliptic case, when 7^ < a(3 and —2 <TrAf < 2. In terms of the complex 
coordinate z = -^{q + ip), the Hamiltonian in reads 

c c 1 1 

H = -z'^ + -i^ + hzz, with 6 = - (a + /?) G R, c = - (a - /5) - 17 G C. (12) 



and A = y|c| —6^. We shall write M(c,b) for the matrix M constructed via fTI Ul -lfT ^ . 
whenever IP' 7^ 



c 
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We will make use of the following convenient decomposition of a general hyperbolic 
matrix M (TrM > 2). We first introduce some notation. For /i G IR+ we define: 

D{^i) =^ M(c=_i^,6=o), 5(/i) =^ M(c=-^,6=o), R{^i) =^ M(c=o,fe=-;.). 
Clearly, -D(/i) is hyperbolic, with the q and p axes as unstable and stable axes. -B(/i) is 
also hyperbolic, with eigenaxes forming angles ip+ = | arg(— ic) = f = — V^- with the 
horizontal. -R(/i), on the other hand, is just a rotation of angle /i and hence elliptic. Any 
hyperbolic matrix M(^c,b) as in (fTUj) can be decomposed as: 

= gD(A)g-\ with Q = R{h)B{h2), (13) 

where 61 G [— |, |], 62 ^ R are defined as follows. We denote by 0i G [— f , f] the angle 
between the q axis and the bisector between the stable and unstable axes of M(^c,b), and 
by 02 £ ]0, ^] the angle between the bisector and the stable axis of M(c,fe) (Figure EI). In 
terms of those, one has: 

This last decomposition has the following interpretation. The general hyperbolic map 
M(^c,b) is obtained from the special case -D(A) (A = a/|cP — 6^ > 0) by a change of coor- 
dinates Q yielding a transformation from the {q,p) frame into the unstable-stable frame. 
The unstable (respectively stable) direction is given by the vectors v+ = Q Cq, V- = Q Cp 
(which are, in general, not normalized). Above, we decomposed Q into the transformation 
5(62) which changes the angle between the stable and unstable axis, and the rotation R{bi) 
which rotates the whole frame (Figure Ej). 

We finally remark, for later purposes, that there exists another decomposition: given 
M G SL(2, R), 3!c G C, /i g] - vr, vr] so that 

M = M(£,o)i?(/i). (15) 
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Boost B(d,) 




Special Map D(X) General Hyperbolic map M 

Figure 3: Decomposition of the general linear hyperbolic map M(^c,b) as in ()13|) . 
2.2 Linear quantum dynamics 

In terms of the usual annihilation and number operators a = '^^('J + ^P)) ^ = 
i (a'I'a + aa"!") , the Weyl (or canonical) quantization of if in (jH)) is defined as the self-adjoint 
operator H on L^(R) given by: 



111 f c c \ 

H = -af + -(3p^ + 7- {qp + pq) = h i -a^ + -a^^ + bnj . 

The quantum evolution operator for time t = 1 which corresponds to M(c,b) is then: 



(16) 



M(c,6) = exp < -1- 



(17) 



The quantization of the matrix —M(^c,b) = M(^c,b)R{'^) can be defined as M(^c,b)P = PM(^c,b) 
where P = —iR{7i) is the parity operator. The unitary operators M(c,b), M(^c,b)P yield a 
projective representation of SL(2, R) (which resembles the metaplectic representation). We 
will in most of the paper omit to indicate the ^-dependence of the operators H and M(^c,b) ■ 
Let V = ViCq + V2ep G and let : R^ — R^ denote the translation on classical phase 
space by v. The corresponding quantum translation operator is defined by: 



exp 



h 



{vip - V2q) . 



These quantum translations satisfy the algebraic identity 



v+v' ) 



(19) 



with S = ^ {v2v[ — Viv'2) = —-^v A f', so they generate an (irreducible) unitary rep- 
resentation of the Heisenberg group. For any matrix M G SL(2,R), one trivially has 
MT^M~^ = Tmv This intertwining persists at the quantum level: 



(20) 
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3 Classical and quantum automorphisms of the torus 



3.1 Classical automorphisms and their invariant manifolds 



Consider the torus T = jl? clS 3i symplectic manifold with the two-form dq A rfp. Then 
any M G SL(2, Z) defines a (discrete) symplectic dynamics on T in the obvious way. We 
are interested in the case where M is hyperbolic: the corresponding dynamical system is 
then an Anosov system |AAj . The stable and unstable manifolds of any point x G T are 
obtained by wrapping the lines with slopes s± passing through x around the torus. We 
present here some properties of these manifolds that we will need in subsequent sections. 

A simple example we will use for numerical illustrations is the so called "Arnold's cat 
map" |AAj 

'21. 

(21) 



M, 



Arnold 



1 1 



Its Lyapounov coefficient is Aq 



log (^) 



0.9624. The stable and unstable manifolds 



of the fixed point x = are depicted in Figure ^ 




Figure 4: The stable and unstable axes through of the map MAmoid wrap around the 
torus at infinity. We have only represented the first six occurrences. 



For any hyperbolic matrix M, the slopes and s_ of the unstable and stable di- 
rections are quadratic irrationals [i.e. the solutions of a quadratic equation with integer 
coefficients). It is well known |Khj that any quadratic irrational s satisfies the following 
diophantine inequality: 



3C(s) > 0, VA; G Z, V/ G N*, 



k 



>Cis)- 



\ls 



k\ > C{s)j. 



This means that quadratic irrationals are poorly approximated by rationals, in the sense 
that, to get an approximation with an error e, you need a rational with a denominator of 
order at least e^^/^. 
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This inequahty will be used in the following manner. Consider the eigenvectors v± of 
M(c,fe) defined as = Qcq , f„ = Qcp (with Q the matrix defined in Eq. (fT^V As usual, 
their dual basis u± (defined as v+ ■ u+ = 1, f+ ■ n_ = 0, etc.) can be used to express the 
coordinates of a point x in the basis v±: 

X = q'{x)v+ + p'{x)v-, with q'{x)=x-u+, p'{x)=x-U-. (22) 

We call d{x,Zl) the distance between a point x G and = \ {0}, and we will 
estimate it for x on the (un) stable axis: 

3C > 0, Vx e Rv±, d(x, Zl) > ^^^^ 
To prove this, first note that, for any G Z s.t. Ug ^ 0, we have 



|p'(?7,)| = \n ■ = \u^^p\ 



Ug + Up 



> £i£±lbfzi^, (24) 
\ng\ 



where we have used the fact that U-^g/u-^p = s+ is a quadratic irrational. Interchanging 
Ug and np, we obtain a first set of inequalities: 

Lemma 1. There is a constant C (depending on M) such that, for any integer lattice point 

C C 
|p'('^)| > "n — u ^'^^ — "n — n"- 

We can now prove ()23|) as follows. For each x G Mw+, there exists an n G Z^ so that 

■ ^ C± 



d(x, Z^) = ||n — x|| > — — > 



Since, obviously, — x II < I/V2, (I2S1) follows easily. 

We will in addition need a slightly refined statement. If the lattice point n 7^ is 
in a sufficiently thin strip around the unstable axis, it satisfies ||p'(n)f_|| < 1/2 < ||?^||/2, 
which implies the lower bound |q''(^)| > 2\\v+\\ • Together with the above lemma, this entails 

\p'{^)\ > \q'{n)\ ^ certain C[. Interchanging p' ^ q', we see that the same inequality 
holds for points in a sufficiently thin strip around the stable axis. Outside the union of these 
strips, this inequality can be violated by at most a finite set of lattice points; therefore, 
upon reducing the constant C[ we obtain the main technical result of this section: 

Lemma 2. There exists a constant Co > (depending on M) such that, for any integer 
points n ^ m of the plane, their coordinates along the (un)stable directions satisfy: 

\q'{n)-q\m)\> ^\ (25) 
\p\n) — p'[m) I 

These inequalities precisely control the sparseness of the lattice points inside a strip 
around the unstable axis: the narrower the strip, the farther successive lattice points have 
to be from each other. 
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3.2 Quantum mechanics on the torus 

We recall as briefly as possible the basic setting for the quantum mechanics of a system 
with T as phase space, as well as the quantization of the automorphism M, referring to 
|HB| IDE| IBouDBj and references therein for further details. In order to define the Hilbert 
space associated to T, we first consider the translation operators Ti = T(i o), T2 = T(^o,i), 
which satisfy T1T2 = e~^/'^T2Ti as a result of p9|) . So for the values of h defined as: 

N = -^eW, (26) 

one has the property Ti, T2 = 0. The Hilbert space L^(M) may then be decomposed as a 
direct integral of the joint eigenspaces of Ti and T2: 

L\R) = pHMfi-^,, 7i^,, = {|^)e5'(M) |fi|^)=e^^^|^),t2|^) = e'^^|V^)}. (27) 

The 'angle' 6 = {61,62) G [0,27rp thus describes the periodicity properties of the wave 
function under translations by an elementary cell. TiNfi is A^-dimensional. 
We can define a projector from iS(M) onto the space TiNfi'- 

Pe= e-'"^'^'-'"^'''^'' TI' t^^ = J2 e"^^'"+'''" t^. (28) 

(ni,n2)eZ2 nGZ2 

The phase (5„ = —nin2NTi comes from the decomposition T„ = e"^"^" T^^^Tg^. 

The Weyl quantization of a function f{x) = ^^^^2 fke^^'"^'^^^^ is an operator on Hn^ 
defined by 

f=Y.hfu,N. (29) 
For \ip) G T^AT^e, its "Wigner function" W^{x) is the distribution implicitly defined via 



{Mm = / f{^) dx, so that W^k) = {^lj\Tk/N\^) (30) 

Jt 

where the W^{k) = J^.e'^^'^^'^^^^ W^(x) dx are the Fourier coefficients of W^. 

Let now M G SL(2, Z), so that A, B, C, D (see Eq. firUjl ) are integers. One then easily 
deduces from fOUj) and that the quantum map M satisfies: 

MPe = Pe'M, with 6' = 6M-^ + 2n^ {CD, AB) . (31) 

The constant shift on the right hand side (RHS) is due to the phases 5„ appearing in 
(j28|) . M will define an endomorphism in 0.^,9 provided 6' = mod 27r, i.e. provided 6 
is a fixed point of the dual map defined in (j3H) . Given a hyperbolic matrix M, such a 
fixed point exists for any |DEj . In particular, for any matrix M the angle 6 = (0,0) 
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(periodic wavefunctions) is a fixed point if N is even, wliile 6 = (vr, vr) (antiperiodic wave 
wavefunctions) is a fixed point for odd. We will always make this choice for our numerical 
examples. 

From now on, we will assume that M = ±M(^c,b) G SL(2, Z) is a fixed hyperbolic matrix 
defining a dynamics on the plane and on the torus. We will therefore no longer indicate its 
dependence on (c, b). We will also assume that h is such that (j^Uj) holds, and for this h we 
select an angle such that 6' = 9. In general, 9 can depend on but we will not indicate 
this dependence. 

4 Coherent states and their evolution 

4.1 Standard and squeezed coherent states 

With the normalized state |0) defined by a|0) = 0, a "standard" coherent state is 

|x)=f,|0), x = (g,p)GM'. (32) 
More generally, we define for each c G C* the "squeezed" coherent states |x, c) by 

|g) = |0,g) =M(5,o)|0), |x,c)=f,|g), (33) 

where the "squeezing operator" M(c,o) is defined by (I17|l . with 6 = 0. Note that, in view of 
(fT^. given M e SL(2, M), 3!g G C, G [0, 2ti[ such that 

M|0) = e^"|c). (34) 

For more details on coherent states, we refer to |Z t IPej . 

To avoid confusion, we will use a tilde for the parameters of the squeezing operator 
M(c,o)5 cind keep untilde notations for the parameters of the dynamics defined by the matrix 

M = ±M(e,6) that are at any rate kept fixed throughout the further discussion. In the 
L^(M) representation, the state \x,c) is a Gaussian wave packet with mean position q. Its 
Fourier transform is centered around the mean momentum p. For any state {ip) G (M), 
we define its Bargmann function as x i— >■ {x, clip), and its Husimi function to be the positive 
function Hg,^ defined on phase space by: 

'Hc,^ (x) = ^^^''^^^^^ , which satisfies / Hc,^ (x) dx = WipWlir^). (35) 

Note that for given lip), the Bargmann and Husimi functions depend on the choice of 
c. Also, the function x ^ {x, clip) is the product of a Gaussian factor with a function 
holomorphic with respect to a c-dependent holomorphic structure. The term Bargmann 
function is usually reserved for the holomorphic factor, but we find it convenient to adopt 
here a slightly different convention. 
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We will need the explicit expression of the (standard) Bargmann and Husimi functions 
of the squeezed coherent state |c): 

/ ni~\ ^ / .#tanh|cn J I ( \\ toa\ 

(a;, 0|c) = , , exp <^ -1 — \ exp <^ -- — — + — — \ . (36) 



Here the unstable-stable frame {q,p) of the symmetric matrix M(^c,o) is easily seen from the 
formulas in Section |21to be obtained from {q,p) by a rotation of angle ip+ (Figure EI), and 
the widths are given by 

Aq^ = - Ap' = - (37) 

(1— tanh|c|) (l + tanh|c|) 





Ap 



Figure 5: Modulus square of the Bargmann function of a squeezed coherent state |c), as 
given in (jHUj) . The inverse Planck's constant = 1/h = 40, and the squeezing parameter 
c = — i|c|e~^''^+ with |c| = 0.962, = 32° (this corresponds to Ci for the map MAmoid)- 
(a) Three dimensional picture, (b) Typical size and orientation of the distribution: ellipse 
"supporting" the distribution. 

Standard and squeezed coherent states on the torus are defined to be the images of the 
previous coherent states by the projector Pg. We use the notation: 



\x,c,6) = Pe\x,c) e HN,e- 
These states are asymptotically normalized: 

{x,c,e\x, c,9) = l + 0{e-^^^^/^) 
and satisfy a resolution of the identity on the Hilbert spaces Ti-N^e |BouDBj : 

dqdp 



(3J 



27ih 



\x, c, 9){x, c, 6*1 



In, 



(39) 
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Similarly as above, one defines for any \ip) G TiNfi its Bargmann "function" x i— >• (x, c, 0\ip) 
(which is actually a section of a suitable line bundle over T, i.e. a quasiperiodic function 
on M^, but this shall not interest us here), and Husimi function 7ic,'4)fi{x) = N \ {x^ c, 0\iIj)\^ , 
a bona fide function on the torus (of which we omit to indicate the A^-dependence) . 



4.2 The evolution of coherent states 

Before turning to quasimodes, we need to study in detail the quantum evolution of the 
squeezed coherent state |c, 6') which is given by |t; g, e) = M*|£, 9),teZ. We will extend 

this notation to any real time, by \t; c, 6) =^ P6»e~'^*''''|c). Due to (jSl|), the states \t; c, 6) are 
again squeezed coherent states (up to a global phase), so this evolution defines a time flow 
c(t) on the family of squeezed coherent states centered at the origin. All squeezed states 
at the origin have even parity: P\c) = |c), so that the evolution of |c) through the map 
MP is the same as through M (yet, these two maps might require different values for 6, 
see Eq. (jSH)). 

It will turn out that \t] c, 6) will be most simply described if the initial squeezed state 
\co,d) at time t = is well chosen in terms of the decomposition (fT^ . Defining, with the 
notations of (fT^ - (fTij) . cq = —b2e~'^'^\ it is easy to check that |co) = e~''^i/^Q|0) since 
M{co,o) = R{hi)]3{h2)R{-hi) and ^(-6i)|0) = Then, with M = QD{\)Q-\ 

M%) = e-'^^/2QD(At)|0), (5o|M*|5o) = (0|D(At)|0) = , \ ^ G M+, (40) 

y cosh(At) 

so the overlap (co|M*|co) is real positive for all times. 

For later purposes we note that, defining, for s eM, Cs E C,as E [0, 27r[ by 



|co)=e'^=|cs), (41) 



(see ()34|) ). it is clear that {cs\e^^^^^^\cs) is real positive for all t. In fact, it can be shown 
that the Cg are the only values of c with this property. Among all s, s = maximizes 
|(0|cs)|^, so |co) is in a sense the most localized state among all \cs). 

In this paper, we will almost exclusively build quasimodes from coherent states with 
"squeezing" cq; this choice is made for pure convenience, and our main semiclassical results 
apply to more general squeezings as well (see Section and Appendix I1U.2|) . 

Before turning to |t; c, 9) G TCN,e, we first describe the evolved state |t; cq) '= e~'^^^^^\co) G 
L^(M), by studying its Husimi function on the plane, as defined in It will be con- 

venient (but again not absolutely necessary for our results, see Section 110. 2|) to adapt the 
choice of c in the definition of this Husimi function to the dynamics M by putting c = cq. 
One then computes 



-H^oA^) = = " o,. ' ^= - (42) 



2 



Co\n\t; Co) m^nQD (At) gtg|o) (0|T^_,^.Z} (At) |0) 



2 



2 



2nh 27ih 2-Kh 
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It is now natural to use the coordinates {q',p') = Q^^ {<l,p) ^ attached to the unstable- 
stable basis {v+^v_) (see Eq. (1^ ). In terms of these, the Husimi function is a Gaussian 
drawn on the unstable and stable axes: 



27r/icosh(At) 



exp 



(1 



12 



/2 



with 



12 



2h 



1 - tanh(At) 



~ he 



2Xt 



Ap 



/2 



Ag'2 Ap'- 



2h 



-2At 



1 + tanh(At) 



Aq' 



/2 t-+oo 



h. 



(43) 



(44) 



The Husimi distribution of the evolved state \t; cq) therefore spreads exponentially (with 
rate A) in the unstable direction of the map, and has a finite transverse width a/^. It "lives" 
in an elliptic region of phase space centered on the origin and of area Aq'Ap' ~ he'^^. Due to 
conservation of the total probability, the height of the distribution decreases exponentially. 
We now turn to \t] cq, 6) = M*|co, 9), t > and its Husimi function 



K 



Co,*,' 



(x) = N\{x,co\Pe\t] £o)P 



It is clear from (PH|) that {x,co\P0\t;co) is obtained by summing (up to some phases) the 
translates of the function {x, co\t; cq) into the different phase space cells of size 1 centered 
on the points of (the cell around will be called the fundamental cell JF). Consequently, 
it follows from (j43j) - (j44j) that this function is non-negligible at a point x ^ T only if x 
lies within a distance ^Jfi from a stretch of length Ag' ^ v^e'^* = e^*^*""""/^^ of the unstable 
manifold through (Figure IH))- Here we introduced the Ehrenfest time as 



T 



def 



log/j-l 

A 



(45) 



Since at time | log /?.|/(2A) = T/2, Ag' reaches the size 1 [i.e. the size of the torus), 
it is clear that for shorter times the Husimi function 'Hco,t,8 lives in an elliptic region of 
shrinking diameter v^e^* around 0. 




t=-6 



t=-3 




q 

t=0 



t=3 



t=6 



Figure 6: Husimi function of the state |t; Cq, 6) for the dynamics ()21|) and = l/(27r^) = 
500. One has T ^ 8.37. 

For times larger than T/2, this Husimi function starts to wrap itself around the torus 
along the unstable axis or, equivalently, the support of some of the translates {x+n, Co\t; cq) 
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start to enter into the fundamental cell. The diophantine properties guarantee that the 
branches of the piece of length Ag' of the unstable manifold passing through the origin are 
roughly at a distance 1/Ag' = e"'^*-*""'"''^'' from each other (Fig. EJ. Consequently, as long 
as Ap' << Q-^i^-^/"^)^ i.e. as long as t < (1 — e)T, the main contribution to {x, co\P0\t; cq) 
and hence to the Husimi function Ti.co,t,e comes from a single term {x + n, co|t; cq) for most 
X G JF. We say there are no interference effects. The regime (1 + e)T/2 < t < {1 — e)T 
was studied in |BonDBl] where it was proven that on that time scale the Husimi function 
equidistributes on the torus. 

For longer times t > {1 + e)T, when the area Ap'Aq' occupied by the support of Hcct 
becomes larger than the area of the torus itself, several terms may contribute equally 
to {x,Co\Pg\t;co). In the next subsection we give a detailed control on the onset of this 
"interference regime" up to time 2T for the Husimi function of \t; cq, 9) evaluated at the 
origin x = 0; we shall show that the interferences remain "small" up to the time 2T. 

As a last remark, we point out that the above discussion is symmetric with respect to 
time reversal. For negative times, Ti-co^tfi spreads along the stable direction, reaches the 
boundary of JF around — T/2, and will interfere with itself for t < —T. 



4.3 Estimating the interference effects 

As explained in the introduction, our crucial technical estimate concerns the autocorrelation 
function for the state |co, 9), given by (cq, 6'|M*|co, 9). More generally, we will need control 
on 

(c„ 9\MXcs. 9) = {cs\PgMXcs) = {c,\MXcs) + Ht, s), (46) 

where we separated the contribution of the term n = (0, 0) (the "plane overlap"), from the 
remaining terms: 

Iit,s) = J2 e"'"'^''" (es|tnM*|5s). (47) 

neZ2 

This remainder represents the interference of the evolved plane coherent state with the 
lattice-translated initial state. We will show that these contributions tend to as ^ oo, 
uniformly for all times |t| < 2(1 — e)T, for any fixed e > 0. 
A trivial upper bound is 

\I{t,s)\ < J2 |(n,£.|e-^«*|c,)| ''^'Mt,s), (48) 

and we shall estimate the RHS. Note that we extended I{t, s) in the natural way to real 
times t. The detailed proofs of the estimates below are given in Appendix llU.ll here we 
limit ourselves to explaining the underlying ideas and to an instructive comparison with a 
numerical example. For simplicity, we will concentrate on the case s = 0. 

We define a time-dependent metric on the plane adapted to the Gaussian in PHj) : 

q'{x) \ ^ If p'{x) 



2 \Aq'{t) J 2 \Ap'{t) 
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The RHS of (jlHj) is simply the sum of this Gaussian of height Ht = (coshAt)^^/^ evaluated 
at all nonzero integer lattice points. The diophantine properties proven in Section 13.1 1 pro- 
vide information on the position of the integer lattice with respect to the ellipse = 1} 

and allow us to prove the following estimates: 

• for relatively short times (meaning \t\ < (1 — e)T), all lattice points n ^ are 
far outside the support of the Gaussian so that \\n\\t is large. In fact, the distance 

reaches its minimum for a single point Uo (more precisely a finite number Af of 
points), with \\no\\f > c ^ >> 1. Note that, here and in the following, we write 
f{h) « g{h) when lim^j^o /(^)/5'(^) = 0. Jo(t,0) is dominated by the contribution 
of this finite set of points, given by A/'-fft exp { — ||no||j }, the contributions of farther 
points being much smaller. The precise bound proven in the appendix reads: 

|t| < T ^ |/(t,0)| < 2v^e-^l*l/2 exp |-Co^| [l + Ce^d*!-^)/^] , (49) 

where the constant Co is the parameter of the diophantine equation ()25|1 . and C can 
be computed explicitly (it depends only on M). 

• For times \t\ > T, a large number of lattice points {N't = Aq'(t)Ap'{t) ~ e'^'-'*'"'^^) 
are contained in the ellipse {i.e. satisfy \\n\\t < 1), and their collective contribution 
dominates the RHS of \I{t)\ < MtHt ~ qK-^MA^) . This is indeed essentially 
what we prove: 

T < |t| ^ |/(t,0)| < ^^e"(-T+N/2) [1 + c'e^(T-|t|)/2] ^ (^q) 

where C can be computed explicitly in terms of M. This upper bound becomes of 
order unity for |t| — 2T. 



From the definition ()46p . we have trivially for any time 

\I{tM < (£o,^|£o,^) + (£o|M*|£o) <l + (e-c(^«)/'^) + ^ 



/cosh(A|t|) 

Combining these estimates (generalized to s 7^ 0), one obtains the following proposition: 

Proposition 1. There exist positive constants C , C , C" such that for all times t G M, 
and for all s in a bounded interval 

\I{t,s)\ < Jo(t,s) <min(c;ie"l*l/2, 1 + v^e'^l^l/^ + CV^"/^) . (51) 

This shows that the interferences remain small until times of order 2T. The existence 
of "short quantum periods" for certain values of h (see the introduction and Section IHI 
implies that l{t,0) is of order 1 at t = P ^ 2T for these values of h. This is further 
illustrated in Figured 
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log ll(t,0)l for l/h=N=9349 




Figure 7: Numerical calculations of log |/(it, 0)|, for the map ()2H). The heuristic upper 
bound F{t) (solid line) is defined in terms of the shifted Ehrenfest time T' = log(A^)/A: 
F{t) = - f - e^(T'-t) _ 1 for < t < T' and F{t) = f (t - 2T') + 0.5 for T <t< 2T . The 
horizontal dashed line at \og{h}/'^) gives the order of magnitude of the plateau for t > T' . 



Figure [7| shows numerical calculations of log|J(t, 0)| for values of Planck's "constant" 
N = 9349 9359 and compares them to F{t), which is essentially given by the upper 
bounds (P^ - lj^Uj) . We observe that, whereas (jl^ is close to optimal, the same is not true 
for for most values of N: there is a "plateau" log \I{t, 0)| ~ log(^^/^) for t > T', where 
T' = log(A^) /A is a shifted Ehrenfest time. This plateau can be explained by assuming that 
the phases which multiply the different terms in I[t, 0) are uncorrelated, like independent 
random phases. For t » T, the RHS of (jTTj) could then be replaced by a sum of many 
(~ A/t) terms with identical moduli Ht but random uncorrelated phases, similar to a 2- 
dimensional random walk. The modulus of the sum [i.e. the length of the random walk) 
has a typical value \I{t, 0)| ~ y/JJtHt ~ h}^'^, independent of time: this is indeed what we 
see numerically. 

However, for the value of = 9349, corresponding to a "short quantum period" P = 
19, as discussed in Section |HJ log|/(t,0)| is close to the upper bounds up to 

time P ~ 2T'. In such exceptional cases — crucial in this paper — there appears strong 
correlations between the phases in the sum /(t,0): the random walk somehow becomes 
"rigid" , which makes its total length of the same order as the sum of individual lengths, 
\I{t, 0)1 ~ Jo(t, 0) ~ MtHf. This rigidity can actually by analyzed directly from the explicit 
expression for the phases |FN2j : one first finds that for these special values of Planck's 
constant N = Nk and t in the interval T < t < 2T, the phases corresponding to the 
relevant ~ Aft lattice points are all close to 2d-th roots of unity, where d = (trM)^ — 4 
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(in the example M = MAmoid and Nk = 9349, the relevant phases are all close to unity). 
Then, the sum of these ~ Mt phases behaves like G{M, Nk)N't for At >> 1, and one can 
check that the prefactor G{M, Nk) (a Gauss sum) is bounded away from zero uniformly 
{e.g. (j(MAmoid, 9349) = 1). This explains the behaviour |/(t, 0)| ~ Jo{t,0). This situation 
drastically differs from the case of a "generic" N, where the relevant phases are more or 
less equidistributed over the circle. 

5 Quasimodes at the origin 

5.1 Continuous time versus discrete time quasimodes 

We are now ready to study the quasimodes and (jHl) "associated" with the periodic 
orbits of the dynamics generated by M, as discussed in the introduction. To alleviate the 
notations, we start with the case where the orbit is simply the fixed point (0, 0) e T. The 
rather straightforward generalization to arbitrary orbits is given in Section [7| Note that 
the Ehrenfest time T = -^^^ is in general not an integer: whenever T or T/2 appears in a 
sum boundary, they should therefore be replaced by the nearest integer. 

It will be convenient to also consider slightly modified quasimodes, for which the initial 
state is not the squeezed coherent state |co, 0) as in but rather the following superpo- 
sition of squeezed coherent states: 

Pe [ dte-'"^'e-^^'\co). (52) 
Jo 

The "continuous time" version of the quasimodes defined in (j2)) then reads: 

l^cont^ def f' iHt|g^^ ^53) 

Jo 

= Pe j dte-''^^e-^^'\co). (54) 

Here we introduced, for any G M, to < ^ ^7 the operator 

ti-i 

AoA.* = E^"'"^M\ (55) 

t=to 

and the equality (jH^ follows from a trivial computation. 

These quasimodes can also be decomposed into 4 parts |$™^*), obtained by integrating 
in t over time intervals of length T/2, then projecting the obtained state in TiNfi- A 
remarkable and useful property (derived from Poisson's formula) is that we can recover the 
"discrete time" quasimodes |$^''^'^) defined in (j21) from the "continuous time" ones: 
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Notice that the state in ()52|) is not 27r-periodic with respect to so that the quasimodes 
1$^°°*) depend on the "quasienergy" G M. 

The main reason for considering continuous time quasimodes is that they are easily 
connected with generahzed eigenstates of the Hamiltonian if, which allows to pointwise 
describe their Husimi densities, a task we turn to in Section IHl 

In the next subsection, we start our study of the above quasimodes. We will use from 
now on the notation |$^) in statements that are valid both for l^^'^'^) and 1$™°*) (and 
similarly for l^j,,/,)). 

5.2 Orthogonality of the states |$j,0)n at fixed 
Proposition 2. (i) The states \^j,(j>), j = 1,2,3,4 and satisfy, as h ^ 

^1$,, ^) = |5i(A, 0) + 0(1), ($^|$^) = 2TS,{X, 0) + 0(1), (56) 

where the smooth function Si{X, (j)) is strictly positive for all (p and 0{1) is uniformly 
bounded in (p. In particular these states do not vanish for small enough h and the normalized 
quasimodes satisfy (jS))- 

(a) Furthermore, for all G M, the become mutually orthogonal in the semi- 

classical limit: for all j ^ k & {1, ... ,4} 

lim n{^j,^\^k,<p)n = 0. (57) 

re— >0 

The limit is uniform for all (p in a bounded interval. 
(Hi) Consequently, for all G M, 

„($0|$j,^)„ 1/2, „($0|$erg,0)n ^ 1/^2 and n{^4,\^\oc,4,)n ^ 1/^2- 

Proof, (i) We first give a detailed proof for the "continuous time" quasimodes. Writing 
k = j — i E {0,1, 2, 3}, a simple computation yields (see (jS}) 

Z_l 

(<l>^7|<l>^7*) = Y.Y. f r^,/g-i(t-t'+s-s'+kT/2)0 ^g^|^^^_^H(t-t'+s-s'+kT/2)|g^^_ 

i=0 t'=0 ^ ° 

Using and PHj) this becomes: 
where 

error < ds ds' {— - \s\) Jo{t + k— + s - s' , s'). 
Jo Jo T 

*~ 2 
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Using the bound (j5T|) . one readily finds that the second term is 0{h~^). 

To estimate the norm of l^™'^*), there remains to compute the integral in ()58|) in the 



case i = j, that is A; = 0: 

fT/2 



T/2 



ds (T/2 - |s|)e-'^'^(co|e-^"^|co) 



T/2 (r/2- |s|)e-^^^ 

ds = — 

.T/2 A/cosh(As) 

|sr(A,0,T/2)-SrnA,0,T/2), 



where the (real) functions S^"^^, are defined as follows: 



^r(A,0,r) 



def 



dt 



-it(f> 



T ^/cosh (At) 



5rXA,0,r) 



def 



Itle"^*'^ 



T ^/cosh (At) 



(59) 



The limits of S'j'^°°*(A, 0, r) as r — > cx) clearly exist. We only give the value for Si , the 
most relevant one for our purposes [BaTITj : 



^r*(A, 



def 



lim ^™'^^(A,0,r) 



AV2^ 



(60) 



For fixed A, this function is maximal for = (with value ~ 5.244/A), and decreases as 
y^^^g-'r|</'l/2A ^Qj. _^ crucial property is the strict positivity of this function, for 

all values A > 0, G M. 

The computation of is similar. 

(ii) We now estimate the overlaps ($i_<^|$j^0) for j ^ i, by estimating the first integral 
of (j3Hj) in the cases 3 > A; > 1: 

r.T/2 



T/2 



ds (T/2 - \s\ 



-i(s+kT/2)(/) 



v/cosh(A(s + fcT/2)) 



4a/2 _ A(k-l)T 

< -^e 4 = c>(;i 



k -l , 
4 



Taking into account the estimate of the error in (j58|) , we see that for any i ^ j, the overlap 
i/^cont|^cont^ is bouudcd by a constant (even by C(/i^/^) for |i — j| = 2). As a result, 



< — . 



/<f,cont|^cont\ — ' 



(61) 



This proves (ii). Part (iii) is now obvious. 

To treat the case of the discrete quasimodes, the integrals over time have to be replaced 
by sums over integers. For instance, the expressions defined in are replaced by 



5r(A, 0, r) ''^ V 

V cosh (At) 
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and similarly for 82- The sum Sf^'^'^{X, 0) = limT-_+oo Sf^'^'^{X, 0, r) is also nonnegative for all 
A > 0, G [— 7r,7r]. Indeed, Poisson's formula induces the identity 

5f^^(A, 0) = 5^ Sr^A, + 2^77). 

fcez 

The norms of the discrete quasimodes therefore satisfy an estimate similar to (jSTj), upon 
replacing Sl°^^ by Sf^^'^. The other estimates are identical as for the continuous version. □ 



5.3 Quasimodes of different quasienergies 

We now compare quasimodes 1$,^) of different quasienergies and show: 
Proposition 3. Let 0o be an arbitrary angle in [0,2tt[, and 

TC 

(j)k = 4>o + j;k, k = l,...,2T. 
The 2T quasimodes |$0j.)n become mutually orthogonal in the semiclassical limit: Wk' 7^ k, 

This is an immediate consequence of the following finer estimate: 

Proposition 4. Let / C M 6e a fixed bounded interval. There exists a constant C > such 
that, given any semiclassically vanishing function 6{h) and n G Z*, if (j),(f)' & I, and if the 
phase shift A(j) = (p' — (p satisfies |A0 — ^| < 6{h)^^, then we have, for small enough h, 

n{^^'\<^^)n<C{9ih) + ^). 

Proof. As before, we write the proof for the continuous time quasimodes. The overlap 
(|<|)Cont|^cont^ givcu by an expression similar to Using the estimate for I{t,s), 

we obtain 

.2T Vcosh(As) A0/2 

where we introduced (f) =^ ^-j^- This integral is bounded above by so that for 

a phase difference bounded away from zero {i.e. |A0| > c > 0), the scalar product of 
the normalized states is ,j>'\^ cj)) n = C(T~^). We are however more interested in the 
case where A0 is ^-dependent and semiclassically small: A0 — > 0. Inserting | sin{A0(T — 
|s|/2)} — sin{A0T}| < |s|A0/2 in the integral and using (j56p . we get for — > 0, A0 — 0: 

/$cont|$r) = sin(TA0) 
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The first term can be as large as 1, for A0 << T ^. It will also be large for values 
A0 = ^(ii±iM with n an integer, \n\ « T, where it takes the value ±1/TA0. At the 
opposite extreme, the term vanishes for A0 = n a nonzero integer, and close to this 
value it behaves like (-1)"^(A0 - □ 

We are now set to analyze, in the next subsections, the phase space distributions of the 
quasimodes \^^)n and of their components. 



5.4 Localization of |$ioc,(/))n near the origin 

Recall that |$ioc,(/)) = |'^2,</-) + I'^^s,.^)- We will show the following: 
Proposition 5. Let G M. Then, for any f G C°°(T2), 

lim „($ioc,<^|/|$ioc,<^)n = /(O) and \im f{x)Hco,\oc,ct>,e{^)dx = f{0). (62) 

where 'Hco,ioc,ci,,e{x) = (x, Cq, 6'|$ioc,</,)nP is the Husimi function of \^\oc,4>)n- H follows 
that the semiclassical measures 'Hco,\oc,ii),e{x)dx and the Wigner distribution converge to the 
delta measure at the origin. All limits are uniform for cf) in a hounded interval. 

Using a more physical terminology, one can say that the quasimodes |$ioc,(/))n strongly 
scar (or localize) on the fixed point G of the map M. 

Proof. As before, we write the proof for given by 

rT/2 

\^Z:^)=Pe / dte-'^'\i-~Co). 

J~T/2 

This is a sum of evolved coherent states for times \t\ < T/2. At this maximal time, the 
length Aq' of the Husimi function of e~'^^*|co) reaches the size of the torus. To control 
the contribution of the nonlocalized states at t ^ T/2, we first select a function Q{h) such 
that in the small- limit 1 << Q{h) << T. We then split pieces: 

IKf^) = Pe [ dte-'^%Co)+Pe [ dt e-^<^*|t; cq) 

where t^, =^ [T/2 — 0(^)/A]. From the proof of Proposition |21 it is clear that 

($'1$') ~ 2nSr\X, 4>) ~ TSr\X, 4>) ~ ($ioc,</,|$ioc,</,) when h ^ 0. (63) 
The norm of the remainder |$") is estimated similarly: 

($"1$") _ ^5r*(A, 0) < CQ{h) = o{T). (64) 
A 
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In the interval \t\ < t^,, the elhpses supporting the states |t; cq) have lengths y/he'^^ < 
Q-Bih) _^ Q Considering the disk Dq centered at the origin and of radius e"®*^'^-'/^, the 
Husimi functions of these states are therefore semiclassically concentrated inside D@. We 
will show below that |$™^0)n is also concentrated inside this disk. 

Using (jH!^ and (jEH), together with the obvious \a + 6p < 2(|ap + one finds 

N\{x,Co,e\^^C'Mdx <^f N{\{xrco,eW)\^ + \{xrc,,e\^")\^)dx 

T\De ^ Jl\De 

< N\{xrc,,eW)\^dx + m<^")]. (65) 
^ \Jt\Dq J 

< (66) 

The last inequality comes from the observation that the Bargmann function {x,Cq,6\^') 
is a sum of Gaussians of widths smaller than e"®*^^-* so simple analysis shows that the 
integral in is O (A^exp(— ce®*^^))) . Consequently, holds and yields the proposition 
provided we choose log log << 9 << logA^. For discrete quasimodes, we only need to 
replace 5*™°* by Sf^^^ in the above estimate. □ 

For later purpose, we notice that the previous proof can be applied to the states 

= j'' e'-^^lt; Co), with - | < ti < < ts < ^. (67) 

These states indeed localize at the origin in the sense of equations (jU^ and The same 
is obviously true for the discrete analogues of these states: \^f^t^) = 'Pti,t2\co, 0). Note that 
|$2,(/.) and I $3,,/,) are of this type. 



5.5 Equidistribution of |$erg,(/))n 
Recalling that \^crg,<f>) = I "^1,0) + |'^4,?i) we have 
Proposition 6. Let G M. Then, for any f G C'^{T'^) 

lim „($erg,<^|/|<l>erg,</,)„ = / f{x)dx = lim / f{x)nco,eTg,<i,,e{^)dx, (68) 
h-^o J J h-^o J J 

where Hco,eTg,(i>,e{x) = (x, cq, 6'|$ioc,<^)nP is the Husimi function 0/ |$erg,0)n- It follows 
that the Husimi measure 'Hevg,(i>{x) dx and the Wigner distribution converge to the Liouville 
measure on the torus. The limits are uniform for (p in a bounded interval. 

The states |$erg,(/.) are said to semiclassically equidistribute on the torus. 

Proof. We will use the algebraic structure of the quantized automorphisms in the proof. 
We will drop the index from the notations. It is clearly enough to show that, for each 
A; G Z^, we have 

lim($erg|rfc/Ar|<l>crg) = 0. 

?!,— >0 
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For that purpose, we write 

($erg|Tfc/Jv|$crg) = ($1 1 Tfc/A. | $l) + ($4 |Tfc/iV | $4) + ($1 |Tfe/JV | $4) + ($4 |Tfe/iV | $l) • (69) 

We first estimate the two diagonal terms of the RHS. Using = e"^"^M^^|$3), |$4) = 
g-i(/)T]y[T|^^^ and the intertwining property (j20|l . we get 

($l|ffc/^|$l) + ($4|Tfc/^|*4) = {<^3\fkJ^3) + ($2|T,.J$2). 

Here fc± = M^^k/N e {Z/Nf are of order 1 (see below), so that we transformed the 
"microscopic" translation by k/N (of order K) into "macroscopic" ones. Each term is 
therefore the overlap between the state |$2) or I'^'s) localized in a small disc Dq centered 
at the origin of the torus (cf. Eqs. ()65I67|) ). and a translated state localized in the disc 

Dq^± =^ DQ + k± centered at the point k± mod 1?. This overlap will consequently be small 
provided k± is sufficiently far away from the integer lattice. We prove this fact using ()22j) : 

k+ = q'{k) e^VN v+ + p'(k) e-'^VN v_ = C(N)q'(k)v+ + 0{h'')Y^, 
fc„ = q'{k) e-^VN v+ + p'(k) e^VN v_ = C(N)p'(k)v_ + 0{h^)Y+, 

where 27re-^ < C(N) < 27re^ since T/2 is the closest integer to |ln/i|/2A. Now, q'{k) ^ 
7^ p'{k) since the slopes of v± are irrational. Consequently, k± are at a finite distance 
from for small enough h, and the disks Dq and -De,± do not intersect each other. We 
can thus estimate the overlap: 

n In 

Jt 

= \ +/ \ n{^3\x,Co,9) {x,Co,9\fkJ<!>3)nNdx. 

lJT\De JDe) 

Using the Cauchy-Schwarz inequality, the first integral is bounded as 

/ n{^3\x,co,6) (x, co,6'|Tfc^|$3)„ N dx 

Jj\D0 



< 




where we used (|65|) applied to |$3). The integral over Dq is treated similarly, exchanging 
the roles of both factors: now the second factor semiclassically converges to zero due to 
the inclusion Dq C (T \ -De,+)- In the end, we get for log log << Q{h) << logN 




„($l|T,/^|$l)„ =r,{^3\M'Tk/NM-' |$3)„ = 0\\^] (70) 



uniformly for in a finite interval. The proof goes through unaltered for the second overlap 
n{^4\Tk/N\^4)n and in fact for any M'^\^ti,t2)n as in (|U7j). leading to: 



25 



Lemma 3. Consider a semiclassically diverging function log|log/i| << O(^) << |log^| 
and /c G . Given a bounded interval, there exists a constant C so that for all in the 
interval 

As a result, the states M'^\^ti,t2)n equidistribute as h 0, which imphes that the 
integral of their Husimi function over a fixed domain of area A converges to A. We now 
use this information to finish the proof of Proposition IHl 

We enlarge the Figure ^ and define the additional state |$5) = e~''^"'"M"'"|$3), which, 
according to Lemma El equidistributes. Now, using the same intertwining property as 
above, we rewrite the nondiagonal terms in the RHS of ()69|) as 



with the vector k' =^ M'^^'^k/N. Each term is the overlap between a state localized near 
the origin {e.g. ($2!) and an equidistributed one {e.g. Tk'l^^))- It is natural to expect 
that they are asymptotically orthogonal. 
To prove this fact, we proceed as above: 



n{^2\fk'\^5)n\ < \ dx Hco,2,e{^) \ dxHcofi^x - k 

V Jj\D{r) V Jj\D{r) 



+ A dxHco,2,e{^)^ \ dx7icQ^5fi{x - k'), (71) 

y JD{r) y JD(r) 

where D{r) is the disc of radius r centered at 0. Using the semiclassical localization of 
|$2)n at the origin and the equidistribution of |$5)n, we find 



limsup I ($2|Tfc/|$5)„| < v^r. 



Since this is true for any r > 0, lim^-^o |„('^'2|^fc'|'I*5)n| = 0. We now control all the terms 
of (j69p and after taking care of the normalizations we obtain Proposition [HI □ 

5.6 Semiclassical properties of l^^) = l^ioc,^) + l^erg,^) 

We now finally consider the "full" quasimode | It is the sum of two states, one localized, 
the second equidistributed. 

Proposition 7. For any G M, (0j holds with r = l,Xo = 0. The limit is uniform for 
belonging to a bounded interval. 

Proof. It is again enough to study „($|T'fc/7v|$)n and to show 

hm „(<l'|ffc/^|<l')„ = ^(1 + 5fc,o)- 
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The results of the previous subsections imply immediately that this reduces to showing 

lim n("^'loc|rfc/iv|'^'erg)n = 0. 

This in turn is proven as in the previous subsection through the use of the Cauchy-Schwarz 
inequality and cutting the integral over T into the integral over a small disc around the 
origin and an integral over the complement (see ()71|) ). □ 

To conclude this section, let us remark that the semiclassical properties of the various 
quasimodes we introduced are not altered if we replace T in the sum or integration bound- 
aries by an integer that differs from it by a finite amount, bounded as h goes to zero. This 
will occasionally be useful in the sequel. 

6 Pointwise description of the quasimodes 

In the last sections, we showed that the Husimi and Wigner functions of the quasimodes 
\^4>)n converge to the measure ^-^^ in the semiclassical limit. The crucial tools of the proof 
were, on the one hand, precise estimates of the overlaps (cq, 6\M^\co, 6) (obtained using the 
diophantine properties of the invariant axes) , on the other hand the algebraic intertwining 
between M and the quantum translations. 

Still, it would be interesting to know the speed at which this convergence takes place, 
or to compute more refined "indicators" of the localization of the quasimodes. 

In this section, we will use a more "direct" yet slightly more cumbersome route which 
will yield more precise information on the phase space distribution of the "continuous time" 
quasimodes. The main step of this route is the pointwise description of the Bargmann and 
Husimi functions of |$™°*). This description will then provide an estimate of the speed of 
convergence to the limit semiclassical measure; at the same time, it will allow us to compute 
alternative localization indicators, like the L'^-norms of the Husimi functions. The pointwise 
estimates will also uncover the "hyperbolic" structure of the Husimi functions near the 
origin, a structure already emphasized by several authors for finite-time quasimodes |KHt 
IWBVBj and for spectral Wigner and Husimi functions |RUdAj . 

6.1 Plane quasimodes 

Our final objective is to estimate the Bargmann function (x, cq, ^l^™"^*) for x E T the 
fundamental domain. For this purpose, we start from quasimodes of the Hamiltonian H: 

\^^,t) = j\se-''^%-'^^'\c,). (72) 

The torus quasimode 1$™"^*) is obtained by projecting |\I^<^,t) onto 71^^ (cf. Eq. (jSll))- In 
this subsection, we will study the Bargmann function of the plane quasimode |^I/0,t)- 
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Using the rescaled variable Z ^= this function is given by the following integral: 

^,^t{x) {x,Co\^,,t) = e-l^l^ r ds -^^e^^^-*-^-. (73) 

J-T V cosh As 

Through the change of variables U = Z^(l — tanhAs), and using the parameter /j, =^ 
l/4 + i0/2A, this integral may be rewritten as 

with the boundaries Uq = ^^(1 -tanhAT) ~ 2^2/12^ Ui = ^^(1 + tanhAT) ^ 2^2(1 - h^). 
This function satisfies the following symmetries (with obvious notations): 

^<^,t(^) = ^0,t(-^) = ^-0,r(i^). (75) 

The hyperbolic Hamiltonian H admits no bound state in //^(IR), but for any real energy 
E = —h(f), it has two independent generalized eigenstates, distinguished by their parity. 
In the limit t — > cx), the quasimode I^E'^,*) converges (in a sense explained below) to the 
even eigenstate, that we denote by {'^Ij^^'^"'''). From the identities H{x) = Xq'p', H = 

XQ^^^^Q~^, the Bargmann function of |v|/^'^*"^"'^^ can be expressed in terms of parabolic 
cylinder functions |NVH iBaHTFj : 



(x,Co|vl/Jr''^")) = C^e-l^l' {D_i+2^(2Z) + D_i+2^(-2Z)}. (76) 

The normalization coefficient = tt (2^' cosh(7r0/A)r(yU + l/2))~^ can be computed from 
the value at Z = 0. For fixed (p and h small, this Bargmann function takes its largest values 
close to the origin (where it takes the value S'™'^*(A, </>)), and is otherwise concentrated along 
the hyperbola {q'p' = —h(j)/X}, which is the classical energy surface {H{x) = —h(f)} (see 
below and Section IfT^ for more details). The Husimi functions of two of these generalized 
eigenstates are displayed in Figure |H1 in terms of the coordinates (Q', P') = ^'^J^r^ - 

From the integral expression (fTSj) . we see that the Bargmann functions of [v]/^'^^'^"'^^ and 
I'^iP^t) are semiclassically close to each other: 

'^(P,t{x) — '^l^^^'^\x) = 0{h}^'^) uniformly with respect to x and 0. (77) 

This equation together with (ffF)|) yields a uniform approximation for '^^^t{x)- One 
cannot simplify this expression in the central region {x = 0{-\/h)}. On the other extreme, 
one can obtain asymptotic expansions for (fTIj) in the region {|a;| >> \/h} {{\Z\ » 1}). We 

will give formulas uniformly valid in the "positive sector" S+ '== {Z \ arg(Z) < ^(1 — e)}, 
where e > is fixed. The symmetries (j75j) then allow to fill the remaining three sectors 
(around the angles 7r/4 + n'7r/2, the function is exponentially small). 
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Figure 8: Husimi functions of two generalized eigenstates l^'^''^"''), in the coordinates 
{Q\ P'). The densities are plotted in linear scale, the contour step depending on the plot. 
The classical energy hyperbolas are drawn in thick curves. 



Expanding the last factor in the integral ()74j] into powers of l/Z"^, we get a sum of 
incomplete Gamma functions jBaHTFl Chap. 9]: 

fj ^f/^e-u (l - = (7(/i,Uo)-7(/i,Ui))+^^^(7(/^+l,Uo)-7(/^+l,Ui))- 

These gamma functions have simple asymptotics in two regimes: 

• for Uq « 1 << Ui, that is, x G iS+, ^/h << \x\ << they yield 

This asymptotics also holds for the Bargmann function of l^^^^^^"^) in the sector 
\Z\ » 1, Z G 5+: it indeed corresponds to known asymptotics of the parabolic 
cylinder functions D_i+2^t [BaHTF, Chap. 8]. This gives for the Husimi function: 



,<p p' 



(79) 



27rh 2AV^ W-¥\ 

For fixed q' » Vh, the p'-Gaussian of width ^/h is centered on the point p' = 
—h(j)/Xq', that is on the classical hyperbola. The function decreases as ^ along the 
"crest" . 
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in the region |x| >> x G 5+, the Bargmann function is "dominated" by the 
coherent state at time T: 



T(a;) = TTTr-nT7 r^e W-^-"^e"'^2^^ — U + C + Of^' 



^0) 



The crossover between the decay and the e~^~ decay is governed by the function 
7(/i, Uq), with Uq ~ hq'^/2 varying from small to large values. 

6.2 Pointwise description of the torus quasimodes 

Using the results in the last section, we will now derive semiclassical estimates for the 
Bargmann function of the torus quasimode 1$™"^*): 

= E^''^''"^^^.T(x + n), (81) 

with the phases i}{x, n) = n ■ 6 + i6n — itcNx A n. 

From now on we restrict x to the fundamental domain JF. We will split the above sum 
between a few "dominant terms" and a "remainder" , which we then bound from above by 
using similar methods as in Appendix llO.ll We will only provide a sketch of the proof. 

From the last subsection, we know that the function "^^^xlx) is concentrated along the 
hyperbola {p' = —h(f)/Xq'}, which is itself a/^-cIosc to the stable and unstable axes. We 
therefore define two strips B^, Bg around these axes: 

Bu = \xeB?, \p\x)\ < 2v^and \q{x)\ < , 5, = {q ^ p'}. 



9VhT 

We call 5 =^ i?„ U Bg the union of these strips, Sq = B^Ci Bg the "central square" and 
i?J and B'^ their periodizations on T or JF. The coefficient Co/9 in the above definition is 
chosen such that i?„ (resp. Bg) does not intersect any of its integer translates (see Eq. ((23) )■ 
As a consequence, for any x E T the intersection between the lattice x + 7? and B^ (resp. 
Bg) is either empty, or it contains a single point noted x + n^^x (resp. noted x + ng^^)-, with 
nu/s,x £ These (possible) points define our "dominant terms" in ()81|) . The remainder 
thus consists in the sum over n El? such that (x + n) ^ B. In order to state the pointwise 
estimate, we define the "modified characteristic functions" Xu{x), Xs{.x) on as 

Xu{x) = e"'('''°"''') if X G -BJ, otherwise 
Xs{x) = e'''('''°=''') if x G -Bj \ Sq, otherwise. 

(this definition is consistent: n^^x is well-defined iff x G -B^). The slight asymmetry between 
Xu and Xs will prevent double counting for x in the central square. 
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Proposition 8. The Bargmann functions of the quasimodes 1$™"^*) have the following 
expression, uniformly for x ^ T and (p in a bounded interval: 

{x, Co, ^1$^°°*) = Xu{x) {x + n„,,, 5o\^^,t) + Xs{x) {x + n,,,, £o|^</,,t) + 0{n"''T^'^). (82) 
On the RHS, |\I^0,t) "^o?/ be replaced by \-qj^^^^'^'>y 

Notice that "^fcp^Tix) at the "edge" of 5„ or B, is of order (^(^i^/^T^/^), so that the above 
estimate of the remainder is sharp. 

This equation gives a precise information for x G -B^, but also a nontrivial upper bound 
in T \ B'^ . It implies that the Bargmann (and Husimi) function of 1$™"^*) is concentrated 
along (a portion of) the periodized classical hyperbola, itself asymptotically close to the 
invariant axes (see Fig. Eland compare with Fig. IH)). These features were not visible in the 
framework of Section El 




Figure 9: Husimi functions of quasimodes |$^'^o) (^^^ -^D linear scale; center: logarithmic 
scale) and |$^'^2Ao) (^igh^> logarithmic scale) of the map MAmoid = 500). 

Sketch of proof. We have to find an upper bound for the sum J^nez'^ x+rv^B l^</',r(a; + '^)|- 
We first consider the points x + n in the sector 5+; since they satisfy |x + r;,| >> the 
Bargmann function is described by formulas (j78|) - (j80p . As in Appendix llO.il we split the 
region B into a union of strips parallel to the unstable axis, of width 5p' = Vh. The 
results of Section jSUl imply that two points (x + n), (x + m) in such a strip are separated 
by at least \q'{n — m)\ > Coh~^^'^. Summing the estimates ()78l80p in these strips, we obtain 
the (x-independent) upper bound 0{-\/hT^/^) for points in iS+. From (f7^. the sum over 
the three other sectors leads to the same bound. □ 

6.3 Controlling the speed of convergence 

Using the pointwise formula ()82|1 . we can now directly compute the Fourier coefficients of 
the Husimi function of l^^""^*): 

^co,*--(fc) = / rfxe2'""^'^N|(x,5o,^|$r')l'' keZl 
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We will prove the following estimate: 

Proposition 9. The Fourier coefficients of the (non-normalized) Husimi function for the 
quasimode 1$^°"^*) satisfy, uniformly for (p in a bounded interval and k (^1? , \k\ < : 

H^,^.o.4k) = Sr'W <P)T (1 + 4,o) + 0{Vt). (83) 

This formula yields at the same time the norm of 1$™"^*), the convergence of the normal- 
ized quasimode to the measure but also the remainder 0{T~^^'^) in this convergence 
(which we could not obtain with previous methods). We do not know whether this esti- 
mate is sharp; in any case, we believe that the remainder cannot be smaller than 0{T~^). 
Using the same methods, we can show that the remainder in the convergence of 
to its limit measure behaves as F{k)T^^, with a function F{k) ^ 0. 

Proof. From Eq. (jH^ . we split Hc^^^^^x) into 3 components: 

^diag(^) = N (|x«(x)*^,T(x + n,,.)P + \Xs{x)^M^ + ns,.)\') (84) 
'M}''"'''{x) = N {xu{x)Ts{x)'^M^ + nu,.)%^ix + n,^,) + c.c. ) (85) 

n~ix) = Oih-'/'T'/') (|x„(x)vi>^,r(x + n„,,)| + |Xs(x)vi><^,T(x + n,,J|) + 0{Vf). 

(86) 

We will show that the integrals over of the "remainder" and the "interference" compo- 
nents are (9(T^/^), while the integral of ^^^^^^^-[^^^^(^x) yields the dominant contribution in 

dHHl). 

The integral of 7-^'''^™^™ on is easy to treat. It involves (ix |\E'0^r(x)|, which we 
estimate by using the asymptotics (|7H|) in the domain x G B, \x\ >> y/h. This yields 
J^dx\^^,Tix)\ = 0{h^/'^T~^/^), so the integral of T^™'! is an olVT). 

Homoclinic intersections To understand the "interference component" 7Y™*°''^(x), we 
have to describe a little bit the set (Sj fl Bj) \ Sq. It is composed of a large number of 
small "squares" surrounding homoclinic intersections (some of them are clearly visible in 
Fig. El). Each of these squares is indexed by a couple of (nonzero and nonequal) integer 
vectors (n„,ns) (finitely many such couples correspond to an actual square in B'^): 

Sqn^,n. = {Bu-n^)n{Bs-ns) = {\q'{x) + q'{ns)\<2^/Kf, \p' {x) + p' {n^)\ < . 

Since we have excluded the central square, one can use the asymptotics (f7S|) for \E'^_r(x -|- 
riu/s)- The integral of |?i'"**^''^(x)| on Sqn^^ns is then smaller than 



C / dq dp 



^y\q'{x + nu)p'{x + ns)\ 
which admits the upper bound 



^/\q'{nu-ns)p'{n,-nu)\ Vk'K-^^)l b'K-^«)| 
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We now want to sum the RHS over all homoclinic squares in B . To compute the sum 
over 1/q' (resp. l/p'), we consider the squares as subsets of Bu (resp. B^), which orders 
them along the strip. Two successive squares do not overlap, so their centers in Bu (resp. 
in Bs) satisfy \6q'\ > A^/hT. As a result, the total number of squares is less than and 
summing their contributions we get 

Notice that we ignored the phases present in 7Y™*°''^(x), as we had done in Section lOl to 
estimate J(t, s). 



Diagonal contribution We now finish the proof by computing the integral 
[ rfx e2''^"^^7^'^'^s(x) = N [ dxe^''^"^'^ |^<^_t(x)P. 

Jb 

The wedge product 2Trx A k is rewritten fc^g' + kgp' in the adapted coordinates. If A; 7^ 0, 
then ku = 27rw+ A k, kg = 27rf_ A k are bounded away from zero (cf. Section IXTj) . 

We give some details for the computation of the integral in the positive sector S+. Let 
Q{h) be a semiclassically increasing function s.t. 1 << Q{h) « T^/^. The integral of 
T^diag ^YiQ central region {\q'\ < Q^/h) admits the obvious upper bound 0{Q^). 

In the region {x G iS+, q' > ©V^}, one can apply the asymptotics (jZHI)- After integrat- 
ing over p', we obtain 

This integral is easy to estimate: 

• for /c = 0, it yields 

• for 7^ 0, it has the asymptotics |BaHTF( Chap. 9] 

^^^^1 log(v^eA:.)| + Oil) = ^n^T + O(log(0fc.)). 
Taking the 3 remaining sectors into account, we obtain the proposition. □ 
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6.4 Husimi function close to the origin and norms 



Besides providing the limit semiclassical measure, the pointwise formula (j82|l allows us 
to compute different indicators of localization for the quasimode namely the 

norms of its Husimi function \Pt\ IN V2] : 



{s>o) \\n^\\s= (^j^[n^{x)Y dx^ 



l/s 



For s = 2, this defines a phase space analogy of the "inverse participation ratio" used in 
condensed-matter physics; in the limit s — > l"*", it yields the Wehrl entropy of the state; 
for s — i> oo, this is sup- norm of the Husimi density. 

Proposition 10. For any fixed oo > s > 1 and (p in a bounded interval, the norms of 
the quasimodes behave in the semiclassical limit as 



C(^,0/A) 

^^"^1 \ogh\ 



By comparison, the L*-norms of a coherent state |c, 6') behave as C'{s,c)h ^+^/* as 
/i — >■ 0, c in a bounded set |NV2] . In the case of the sup-norm, we have a more precise 
statement (see Fig. IS)) : 

Proposition 11. For small enough h, the maximum of l-L^^ \q,<:ont-^^[x) is at the origin for 

I0I/A < 0.5, and C(oo,0/A) = ^'"^"^^s/a!!^/^^''^ • Conversely, for |0|/A >> 1, the maxi- 
mum is close to the points Q' = —P' = ±a/0/A on the hyperbola, and C(oo,0/A) ~ 

Sketch of proof. For any s > 1, the decrease ~ of the Husimi function along the 

hyperbola implies that most of the weight in the integral j^Tii^ont is supported near the 

origin, so that this integral is close to f^i'H^ieven). This yields the proposition, with the 

coefficients C{s, 0/A) given as integrals of parabolic cylinder functions. The statements on 
the maxima derive from known results about parabolic cylinder functions. □ 



6.5 Odd-parity quasimodes 

The connection (|H^ between torus quasimodes |$^°°*) and generalized eigenstates |vl/^'^*"^")) 
hints at a property we have not used much, namely parity. We have already mentioned 
that for each energy E = —h(j), H admits two independent generalized eigenstates, l^^*^*"^"^) 

of even parity, and a second one of odd parity, which we denote by | \E'^°'^'^'' ) . On the one 
hand, the Bargmann function of the latter can be expressed similarly as in Eq. ()76|): 

{x, £o|^r^) = C;e-I^l' {D_i+2,(2Z) - D_i+2,(-2Z)} . 
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On the other hand, as we did for we can build this odd eigenstate by propagating 

an "odd" coherent state at the origin, i.e. replacing the initial \cq) in Eq. (j72j) by the first 
excited squeezed state 

The Bargmann function of the corresponding quasimode |^E'0,r)i is given by an integral 
similar to ()73|). with the integrand multiplied by the factor : this is therefore an 

odd function of x, semiclassically close to to (a;, Co|\E'^°'^'^'') 



(j) / X =0 odd 




Figure 10: Husimi functions of the odd eigenstate I^E^JJLq ) (linear scale) and the torus 
quasimode I^^^^q') (logarithmic scale) for N = 500. Notice the zero at the origin. 

Projecting this plane odd quasimode to the torus through Pg, one obtains a quasimode 
1$^"^'^^) of M with quasiangle (f). Provided one has selected periodicity conditions 9 = 

(0, 0) mod 71, parity is conserved by Pg, so that the Bargmann function {x, Cq, 9\^(p) (resp. 
{X, Co, is an even (resp. an odd) function of x. As a result, these two quasimodes 

are mutually orthogonal. The Bargmann and Husimi functions of 1$^°'^'^'') can be described 
as precisely as for its even counterpart, in particular its normalized Husimi and Wigner 
functions converge as well to the measure with a remainder (9(T~^/^). 



6.6 On the "robustness" of continuous quasimodes 

We want to show that the continuous quasimodes are "stable" with respect 

to a change of the initial state (|co) and |co)i, respectively). One can indeed obtain an 
even quasimode very close to l^*™"^*) by propagating a different initial state \ipo): this 
state needs to be of even parity, sufficiently localized {e.g. a finite combination of excited 
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coherent states), and taken away from a subspace of "bad" initial states. These remarks 
will be made more quantitative in Appendix 110.21 which treats the case where \ipo) is a 
squeezed coherent state of arbitrary squeezing. 

To explain this "robustness", we notice that the operator 



oo 



projects L^(M) onto the 2-dimensional space spanned by |v|/^'^'"^")) and |\l/^°'^'^)). Any even 
state \ipo) G L^(]R) will thus be projected onto Cfp{ilJo)\'^^^^^"^^), with the prefactor 



This prefactor vanishes iff there exists a state {ipo) G L^(M) such that {tpo) = {H + h(f))\ipQ)] 
such \iPq) form a "bad" subspace of codimension 1 inside the space of even states. 

If \ipo) is localized inside a disk of radius Ch}^'^~^^ at the origin, one can describe the 
plane quasimode T'^'r/r.c/ilV'o) as in (f77|) : 

(x, go|^-TX<Al^o) = C^i^o){x, £o|^J""^) + 0{h'/^-') uniformly in x G R\ (87) 

If C(i,{ipo) is of order unity, this estimate shows that 'P^i'rxt/.lV'o) resembles the quasimode 
|^(/),t)- One can then show (as in Section that the torus state PeV'^'^^j.^^lipQ) is close 
to the quasimode C</,('?/'o)|$™''*). 

As an example, consider the case = 0: one can start from any (finitely) excited 
coherent state of the form |co)4n oc M(co,o) ('2^)''"|0) to obtain a quasimode asymptotically 
close to 1$™"^'). On the opposite, the states |co)4„+2 are "bad" initial states, because they 
are in the range of H. 

This discussion straightforwardly transposes to the construction of the odd quasimodes 
l^odd^ starting from odd locahzed states. 



7 Quasimodes on a general periodic orbit 

We have so far described the construction of quasimodes localized on the fixed point 
of the classical map M. We will now generalize this construction to a general periodic 
orbit of M. The associated Husimi densities will be shown to be (semiclassically) partly 
localized on the orbit and partly equidistributed. The proofs require some minor changes 
with respect to the previous case, but no fundamentally new ingredients. 

We consider a fixed periodic orbit V = {xe G ^}J=o '^^ (primitive) period r, in other 
words, for < £ < r, Mxe = xi+i mod 7? and x-r = xq. Note that M'^xi = xe mod Z^, so 

that all Xe, when viewed as points on the torus, are fixed points of M' =^ M'^. Furthermore, 
for all < £ < r, there exist G so that 

Xe = M^xq — me. 
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We will first introduce the discrete time quasimode defined in Q and will consider its 
continuous time analog below: 



T-l 



t=~T 



Letting T be the integer multiple of r that is closest to | ln^|/A, and setting T' = T/t, 
a simple computation yields 



r-l / T'-l 



= ^e-'^^|$f") where |$f = mM ^ e-''^"'^M' Pet,j5o). (88) 

i=o \k=-T' / 

It is easy to see that 

M'Pet, = Pef^,+m,M' = e'^^P^t^.M^ (89) 

where Si = 9 ■ mi + i6mt + '^T^Nmi A xi (see flT^V This phase can partly be interpreted in 
terms of the action along the classical orbit; however, the 6'-term is non-classical, akin to 
the quantum phase due to a pointwise magnetic fiux tube on a charged particle (Aharonov- 
Bohm effect) jEM]- Hence 

|$f ^) = e'^^ [ e-'^'^^Ivi'M Pet,,M^|5o). (90) 

\k=-T' / 

This suggests that l^^'^'^) is a quasimode of quasiangle 0r for M', associated to the fixed 
point of M' . This is basically the content of Proposition^) There is another instructive 
way of rewriting |$^''^'^) which corroborates this idea. For that purpose, we first draw from 



Eq. 

M^'^Pgf^, = (^^^-PeT^^W^ and W^+^PgT^^ = e^('^^-+^^)P0t^,M^'^+^. (91) 
Using this, one can write 

= e'^'T^^Pg^ ( J2 e-'('^"-^-)'^M'M M^|5o), (92) 

\k=-T' / 

where we used PgTx^ = Tx^Pq^ with Og = 9 + 27tN{p£, —qe). A simple computation shows 
that, because X£ is a fixed point for M^, 9£ is a fixed point for the map 9^9' defined in 
(jnH), with M replaced by M'. Consequently, l^^'^'^) is the xi translate of a quasimode for 
M' at the origin with quasiangle (pr — Sr, of the type studied in the previous sections. 

To build continuous time quasimodes, we replace in all the above formulas |co) by 

- r dte-'''^e-^^'\co), (93) 
T- Jo 
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where the "quasienergy" G M is chosen so that 

r0 = r0 - Sr mod 27r. (94) 

Whereas the quasiangle is defined modulo 2tt, the quasienergy is chosen in M. The 
continuous quasimode reads: 

l^cont^ _ gis.r^^^p, 1 /" dte-'^V^^*M^|co). (95) 

' T J-T 

All the above quasimodes can of course in obvious ways be split into a localized and an 
equidistributing part, as before. For both the discrete and continuous time quasimodes we 
have the following estimates: 

Proposition 12. For all <£'<£< t - 1, for all f G C^{T^), for all k e Z, 

= 2T'Si{(PT-Sr,r\) + 0{l) (96) 

hm ($,|/>,)„ = i/(x,) + i/ fix)dx. (97) 
hm „($,,|ffc/^|$,)„ = (98) 

The quasimodes 1$^) satisfy (0), the limit being uniform for (p , cf) in a hounded interval. 

Starting from a pointwise analysis of the continuous time quasimode can 
be performed as well, along the lines of Section 16.21 One should notice that the Husimi 
function of in the a/^- vicinity of a periodic point xi is dominated by the contribution 

of 1$^°^^*); it is concentrated on a hyperbola which depends on the quasienergy rather 
than on the quasiangle 0. 

Proof of the proposition. We write the proof for the discrete time quasimodes only. (jU^ 
immediately implies and (PTjl as a consequence of the results of Sectional To prove 
dHHl) when /c = 0, i.e. the asymptotic orthogonahty of the we write, using (|HH|) and 
(EH) 

T'-l T'-l 

($disc|^disc^ ^ ^ ^ e-'(*^-s^)('^-'^')+'S^-^'(5o|t_.„P,t.^_^,M(^-^')+^('^->^')|5o), 

k'=-T' k=-T' 

SO that 

T'-l T'-l 

|($f=|$f-)|< E i(^oi^-oT^™t,_,M('-'')+^('=-'=')|5o)i 

fc'=-T' k=-T' mGZ2 , s 

T'-l T'-l ^ ' 

k'=-T' k=-T' 

where r = xq — Xi-i', and where we used the estimate Jr{t,0) < C/le''^'*'/^ extracted from 
Appendix 110.11 To prove ()98|) when k ^ 0, one repeats the arguments of Section El we 
omit the details. For continuous quasimodes, the proofs are analogous, using this time the 
same estimate on Jr{t, s). The proof of ((Tj) follows immediately. □ 
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Convex combinations of limit measures We can further enlarge the set of semiclas- 
sical limit measures by taking finite convex combinations of the previous ones. Consider a 
finite set of periodic orbits {Vi, . . . , Vf}, and complex coefficients {ai, . . . , a/} satisfying 
Sf=i Ic^iP = 1- Let \^v^,<j>) be quasimodes (discrete or continuous time ) associated to 
Vi, as defined above, with the same quasiangle (p. We can then combine them into the 
quasimode 

/ 

1=1 

One readily shows along the lines of the proof of Proposition^] that for i j, and for all 

k eT?, one has 

lim n{^V^fk/N\^V,)n = 0. 

This together with (jZj) shows that the Husimi and Wigner functions of converge to 
the limit measure | (^dx + J2i=i l^tiP '^n j • 

8 Scarred eigenstates for quantum cat maps of short 
quantum periods 



We will now slightly extend an argument from |BonDBl| in order to show that the quasi- 



modes we have built and studied in the previous sections are exact eigenstates of the 
quantum map M for certain special values of h and we will prove Theorem ^ 

For that purpose, we first recall a few facts about quantum cat maps [HB^. For a given 
value of = (27r^)~^, every quantum map M has a "quantum period" P{N) defined to 
be the smallest nonnegative integer such that 

MP(^) = e'^^^^^in^^^ for a certain ip{N) G [-tt, 7r[. (100) 

It follows that, if (f) is of the type (p = (pj = "^^^^-^^ , then -pjj^Vti,ti+p{N),(t>j is independent 

of ti, and is the spectral projector onto the eigenspace of M inside TiNfi associated to the 
eigenvalue e^'^j (the normalization factor 1/P{N) ensures that it is indeed a projector). All 
eigenvalues of M on T-Cn,9 are necessarily of that form. 

The function P{N) depends on in an erratic way, and no closed formula exists for it 
|Kej . It satisfies the general bounds 

3C>0, WNeW, hogN -C < P{N) <C NloglogN. (101) 

A 

It is moreover known that, for "almost all" integers, P{N) > \/N |KR,2j . We will now 
give an elementary argument to show that, given any hyperbolic matrix in SL(2,Z), there 
exists an infinite sequence of integers for which the quantum period is very short in the 
sense that it saturates the above lower bound: 

P{Nk) = 2^^^ + 0{l) = 2Tfc + 0(1), (102) 
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where the Ehrenfest time T was defined in 

Let us first recall that, for all A; G N*, one has 

M'' = pkM - pk-i, where pk = _ , Po = 0. 

It was proven in |BonDBl| that, for all k > 1, the integer = GCD{pk,Pk-i + 1) satisfies 

jlogNk = k + 0{l), (103) 

and that 

M'' = I + NkMk, with Mfc an integer matrix. (104) 

We now set Nk — A^^ if A^^ is odd, A^^ — Nk/2 if Nk is even. Choosing the periodicity angle 
6 = (0, 0) when Nk is even and 6 = {it, tt) when Nk is odd (which makes sense, cf. the end 
of Section EH)) , we prove below the following lemma: 



Lemma 4. With Nk, given as above, = e"^i-^j^^ ^ for a certain (p G [— vr, 7r[. 

This means that the quantum period P{Nk) of M on 'HNkfi divides k. Comparing p01|) 
with ()103|) entails that for k large enough, P{Nk) = k and ()102|) holds. 

Proof of the lemma. The case Nk = 2Nk, 9 = (0, 0) was treated in |HBj . We give a different 
proof, which works for both cases. 

From Schur's Lemma and the irreducibility of the Tn/N^i it suffices to show that 

M^,Tn/Nk = on HNk,e, for all n G Z^. Setting Nk = eNk, 6 = e'(7r,7r) and using 
the definition of Pe, Eqs. (fT^ and ((20)), one readily computes 



M'f^/^^M-'Pe = e^-^^'^^^^'^^t./N.TM.nPe 

_ / 1 Ne(nAMfcn)+ee'[(Mfen)i + (Mfen)2+(A'4n)i(Mfcn)2]fi , p 

This phase is trivial if e = 2. In the case e = e' = 1 (that is, Nk odd), one must consider 
the 6 possible values of M modulo 2: in all cases, the phase is trivial. □ 

If we now consider such a value Nk together with an admissible eigenangle 0^^., the 
eigenstates 

P(Affe)/2-l 

\^k)= Yl e-'-^^^^M^Ixcgo) 

t=-P{Nk)/2 

are (discrete time) quasimodes of the quantum map as studied in the previous sections. 
Indeed, as discussed at the end of Section since T and P{Nk)/2 differ by a bounded 
number of terms in the semiclassical limit, we can replace one by the other in without 
affecting any of the semiclassical properties of the quasimodes. One can similarly construct 
eigenfunctions that are continuous time quasimodes. 
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Proof of Theorem^ The previous arguments settle the case /3 = 1/2. To treat the general 
case, we recall that the Schnirelman theorem implies the existence of a sequence of eigen- 
functions \ipk)n of M on 'HN^fi (with corresponding eigenvalues {(t)jk)km) that equidistribute 
as — > cxD. We then construct, for < a < 1: 

\'lljk) = a\^k)n + Vl - a'^\^k)n 

If we show that, for all n G Z^, 

hm n{^k\fn/Nj^k)n = 0, 

a simple computation implies that the \4'k)n satisfy ^ with /3 = We have 

lim n{'Pk\Tn/Nk\^k)n = 1™ f n{'Pk\Tn/Nk\^k,CTg) n + nifk\Tn/Nk\^k,\oc)n] ■ 

The second limit vanishes with an argument as in (fTTjl . whereas for the first, we use the 
further decomposition |<l>fc,erg) = |$fc,i) + l^k^) with \^k,i) = e^'^jkT/2M-T/2|$^ 2), |$fc,4) = 
e~i'?^jk'^/2M'^/2|$]^ 3^ (see ©)• Now, since \^jf,)n is an eigenfunction, we have 

\n{^k\fn/Nj^kA)n\ = \ n{^k\M~''^^fn/N,M'^^''\^k,3)n\■ 
As in the proof of Proposition and more specifically Eq. ()7H) . this tends to with h. □ 

For matrices M of "checkerboard structure", the results of |KRHlMej imply that, given 
an arbitrary sequence of eigenvalues {(pjjkm, there exists a corresponding sequence of 
eigenvectors \ipk) G Ti-N^fi that semiclassically equidistribute. One can then construct for 
the same eigenvalues eigenstates \ipk) satisfying Eq. (jHl). 

The P{Nk) eigenstates with distinct eigenvalues constructed above are of course exactly 
orthogonal to each other, and not just asymptotically as proven in SectionEISl On the other 
hand, two continuous time eigenstates of identical eigenangle (pj but different quasienergies 
4> — (f)' = stt/T, s 7^ become orthogonal in the semiclassical limit. This is also the case for 
two eigenstates with the same eigenangle supported on different periodic orbits V ^ V. 

9 Conclusion 

In this article we have constructed and analyzed a certain class of "quasimodes" of hyper- 
bolic quantized torus isomorphisms, which for certain values of h become exact eigenstates. 
The characteristic property of these quasimodes is that their "quantum limit" , that is the 
weak limit of their Husimi densities, does not yield the Liouville measure, but contains a 
singular component supported on a (finite union of) periodic orbit (s). In our case, this 
singular component has a relative weight (3 < 1/2, less than or equal to the weight of the 
Liouville part. As explained in the introduction, no limit measure of eigenstates can have 
a "larger" singular component. We further conjecture that no sequence of quasimodes {i.e. 
images of the operators 'P-t,t,<I)) can have a more singular limit measure either. 
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The strong scarring of eigenstates exhibited in this paper is directly hnked to the very 
large degeneracies of the eigenvalues of M for certain special values of Planck's constant. 
Therefore, such sequences of eigenstates are very probably absent as soon as one considers 
nonlinear perturbations of the dynamics, for instance = e~^^^^'''^M, for any periodic 
Hamiltonian Hi{x) and e > small enough. Such a perturbation of the classical map 
is known to conserve the uniform hyperbolicity, but destroys the "action degeneracies" 
characteristic of the (linear) cat map. As a consequence the spectrum of the perturbed map 
exhibits Random Matrix statistics, in particular "repulsion" between eigenangles |KMj . 
which forbids degeneracies. 

The precise characterization of some weaker form of scarring for individual eigenstates 
that would remain valid for remains therefore an open problem. Nevertheless, it might 
be interesting to study the phase space distribution of the "nonlinear" quasimodes of the 
type e~"'^*M*|xo, c), for xq a periodic point of M^, which may not be as simple to 

describe as for the linear map. 
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10 Appendices 

10.1 Estimate of the interference term I{t,s) 

In this appendix we prove Proposition ^ For the purpose of Section Q we will at the same 
time give a bound for the more general overlap (t, s G M) 



(105) 



where r ^ T (the fundamental domain) belongs to the lattice (-^Z)^, with G N* and 
where d G [0, 27r[x [0, 27r[ is arbitrary (in other words, Q need not be equal to the fixed point 
of the map (|^. We define 



def 



^ (r + n,g,|e ^"*|cs) 



(106) 



We first consider the case s = 0, t > 0. Since ^Jfi < Ap' < \/2h for all positive times, 
only the points r + n near the unstable axis can significantly contribute. Therefore, we 
subdivide the plane into strips parallel to this axis: the "outer" strips 



V/ > 1, S±i = {x\ai< ±p{x) < ai+i] 
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with ai = Wq + (/ — l)W, and the central strip 5*0 = {x 7^ | p'{x)\ < Wq}- The widths 
Wq, W will be explicitly set below. 

We start by estimating the contribution of the points r + n^Si with Z > 1. Due to the 
diophantine condition ()25|1 . as long as W is small enough, two points in this strip satisfy 
the property |g'(r + n) — q'{r + m)\ > Co/W. Ordering these points according to their 
abscissas: q'{r + rij) < q'{r + rij^i) < q'{r + nj^2), we have for any a > : 



exp |— «g'(r + n^)^} < exp < —a ( 
j&z jez [ ^ 



jCo 

w 



(107) 



The sum on the RHS is a one-dimensional theta function, which has the upper bound 
(optimal for < a small enough): 



(lOJ 



As a result, using PHj) it becomes clear that the contribution to Jr{t) of the points r+n G Si 
is bounded above by 



E 



1 



V cosh Xt 



exp 



p'^r + njY q'{r + nj)" 



Ap' 



/2 



Ag'2 

< y2e-^*/2 e 2^p'- 



1 + V27T 



WAc( 



. (109) 



The estimate pOSp can then be applied to the sum over the strips Si, I ^ 0, to obtain 
(remind |t; Cq) = e~'^*/^|co)) 





J2 \{r + n,Co\t;do)\<V2e-^'/^e'^- 

1^0 r+neSi 



2 + 



/2^Ap' 



1 + 



27rWAq' 



For each time t, we can minimize the RHS with respect to W by taking W"^ = Cq-^^ 
2 



^e which leads to the bound 



Y 5Z K^ + ^,co|t;go)| < 2V2e"^e-^*/2 

Z^O r+ne5; 



1 . J^Ap'e-/' 



:iio) 



Notice that this upper bound is independent of the point r. 

We now estimate the contribution of the strip 5*0, which requires more care, and will 
depend on r. For any point r' 7^ on the lattice sufficiently close to the unstable 

axis, the diophantine property implies |p'(r')| > -fp^fj^y As a consequence, the 
quadratic form appearing in may be bounded inside 5*0 by 



q'{r + ny ^ p'{r + n)^ ^ q'{r + n)^ 



def 



2Ag' 



,12 



2Ap 



J2 



2Ag'2 2D^Ap'^ q'{r + n) 



f,{q'{r + n)). (Ill) 
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The function ft satisfies the scahng property ft{q) = ^r^n f (^^^75^)' '^ith /(g) =^ 

+ q^^- This function /(g) is bounded below for all positive g by the parabola ^'(g) = 
2 + (g — 1)^, so after rescaling we get 



fi „ — At p— 2At , 

V, > 0, /,(,) > S ^ + Ia^S (« - Vae-'VD 

We consider the contributions of the points r + n in S'o such that g'(r + n) > (the points 
with negative g' can be treated identically). We order these points as < go < g^ < . . .: 
each contribution is bounded above by the quantity (cosh At) ^^/^e"^'^^]-*, which is maximal 
for the q'j close to ^/C^e^^/^/D. The diophantine inequality |g^- — g^-^^| > Co/Wq together 
with the estimates p07ll08|) then yield 

This contribution now depends on the rational point r through its denominator D: the 
upper bound increases with D. The full sum Jrit) is bounded above by the sum of the 
RHS in ()ll()j) - ()112|l . For each time t, we adjust the value of Wq to minimize that sum. We 
do not search the exact minimum, but only its order of magnitude. We have to distinguish 
two time intervals: 

• for short times {t « T), the behaviour of ()112|1 is governed by the first exponential 
(since Ap'^ < 2h). We take Wq such that the first exponential in (jllOj) is much 
smaller than that factor, for instance by taking Wq = 2^/C^e~^^^'^ /D. Being careful 
for times around t < T, we find 

< t < T ^ Mt,0) < 2V2 exp [-^^] e-*/^ [l + Ce^^-^^/^] , 

where the constant C is independent on the denominator D. One may replace Ap''^{t) 
by its maximum 2h for positive times. The RHS increases with the denominator D. 

• for times t > T, the RHS of ()112|) is now governed by the factor between brackets, 
and we want to make sure that (IllUj) is not larger than it. Still taking Wq = e~'^*/^ 
leads to the estimate: 

T<t<2T^ J,(t, 0) < ^^ne^*/2 ^1 ^ c'e^(T-*)/2] . 

The constant C is independent of r, so this bound applies uniformly to any point 
a; e T: it yields a L°°-bound for the Bargmann (or the Husimi) function of M*|co, 9). 

The same bounds apply as well to Jrit, s) with s 7^ 0. Indeed, replacing the initial squeezing 
Co by its s-evolved value amounts to dilating the coordinates of the points as g'(r + n) ^ 
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e'*'*^q'(r + n), p'{r + n) ^-^ e~^*ip'(r + n). One easily checks that this dilation does not 
modify the above bounds. 

The negative times are treated thanks to the identity Jrit,s) = J_M-*r(— ^5 -s), and 
noticing that the above bounds only depend on the denominator D, common to r and 
M"*r. 

10.2 Changing the initial squeezing 

We chose from the beginning to construct quasimodes starting from the coherent state |co) 
defined in Section 14.21 The definition was motivated by the positivity property ()4()|1 of the 
overlap {cq\M^\cq), and by choosing the "smallest" parameter c sharing this property. The 
simple expression fl40|) was then used to control the "interferences" J(t, s) (cf. Appendix 
110. and to obtain from there the asymptotic norm of the quasimode f Section 15. 2j) . a 
crucial step for further estimates. Similarly, we also chose to analyze the quasimodes using 
the CQ-Bargmann representation, because of the relatively simple formulas for {x, Co|M*|co) 
(see (gSl)). 

We want to stress (as we did towards in Section lU^ for the continuous quasimodes) that 
both these choices were made purely for convenience, and are not crucial for the results 
of this paper. The construction of quasimodes can be extended in many ways. In this 
appendix, we will consider discrete or continuous quasimodes starting from a squeezed 
state |ci), with an arbitrary (possibly /i-dependent) squeezing ci. We also want to analyze 
these quasimodes using the Bargmann function (x, C2,^|$) for some £2 G C which could 
depend on h as well. 

Proposition 13. The convergence (|7j) holds the above quasimodes, as long as c\ and £2 
stay in a fixed compact set K G C for all h. 

Sketch of proof . For an initial state |ci), the overlap (x, Ci|e~'*^/^|ci), crucial in the cal- 
culation of I{t,s), is still given by closed formulas. We only give it for the simpler case 
x = 0: 

(ci|e-™/^|5i) = (5;|D(At)|5;) = ( cosh(At) + iR(c;) sinh(At))"^/^ 

where \c'i) cx Q~^\ci) and R{c) = — 3?(c)^i^^|||p^. In general, this overlap is therefore not 
real. However, it still decreases exponentially fast with time, and its average 

5i(gi,A,0)= / dte-''^*(ci|e-^*^/'^|ci) 
Jr 

can be easily related with S'i(A,(/)) through a change of variable. One gets S'i(ci, A, </>) = 
a/ cos(Ari) Si(A,0) with the 'complex time' ti = arctan{i?(c']^)}/A. 
For X 7^ 0, the expression for (x, Ci|M*|ci) is more cumbersome than ()43p . Yet, it is still 
a Gaussian having an elliptic profile of width ~ \/h, length ~ v^e^'*' and height ~ e~'^'*'/^, 
and its long axis is asymptotically lined up with the unstable direction for t 00. As a 
result, the results of Sections l4.3H5.HI still hold (replacing S'i(A,0) by Si{ci, \,phi)). The 
localization property ()65p holds as well, even if one replaces in the bras cq by £2, as long 
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as C2 remains bounded. The rest of the proof to obtain ((Tj) f Sections I5.5H5.6|) goes through 
unaltered. □ 



Following the Section Ifj.fjf the plane quasimode "P^T^T.^ki) '^^^ t)e analyzed point- 
wise through the estimate (|H7j): one now has explicitly C(^(|ci)) = e~'^^^ cos(Ari). One 
may replace Cq by C2 in that estimate. As opposed to Eq. (fTBj). the Bargmann function 
(x, C2\^'^^"'^'^^) is not given in terms of cylinder parabolic functions. Yet, its behaviour "far" 
from the origin will be similar to ()78|) . As a consequence, the pointwise estimate (jH^ (with 
Co C2 in the bras) will apply to the torus quasimode -Pe'^^l°r*T,0|ci) as well, upon taking 
the prefactor C^dci)) into account and replacing in the bras cq C2 on both sides. The 
estimates of Sections l(i.3Hfi.4l may be generalized as well to the present case. 
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